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The use of high quality hexagonal boron nitride (hBN) as a dielectric material has made possible 
the realization of graphene devices with very high mobility. In addition hBN can be made as thin as 
few atomic layers and, as recently demonstrated experimentally, can be used to isolate electrically 
two graphene layers only few nanometers apart. The combined use of graphene and hBN has 
therefore opened the possibility to create novel electronic structures. In this work we study the 
"hybrid" heterostructure formed by one sheet of single layer graphene (SLG) and one sheet of bilayer 
graphene (BLG) separated by a thin film of hBN. In general it is expected that interlayer interactions 
can drive the system to a spontaneously broken symmetry state characterized by interlayer phase 
coherence. The peculiarity of the SLG-BLG heterostructure is that the electrons in the layers (SLG 
and BLG) have different chiralities. We find that the difference of chirality between electrons in the 
two layers causes the spontaneously broken symmetry state to be N-fold degenerate. Moreover, we 
find that some of the degenerate states are chiral superfluid states, topologically distinct from the 
usual layer-ferromagnetism. The chiral nature of the ground state opens the possibility to realize 
protected midgap states. The N-fold degeneracy of the ground state makes the physics of SLG-BLG 
hybrid systems analogous to the physics of 3 He, in particular given the recent discovery of chiral 
superfluid states in this system [J. 



Graphene 0| and bilayer graphene H are systems with 
truly unique electronic properties They are ideal 

2D electronic systems in which the conduction and va- 
lence bands touch at single points, charge neutrality 
points (CNPs), at the corners of the Brillouin zone (BZ). 
Around these points the low energy electronic states are 
well described as massless Dirac fermions in SLG and as 
massive chiral fermions in BLG. Recently the use of hBN 
as a substrate has allowed the realization of SLG and 
BLG devices with extremely high mobilities @. More- 
over, hBN films only few atoms thick can be fabricated 
and used to realize graphene heterostructures 0, H| in 
which graphene layers are only few nanometers apart but 
electrically isolated [§-[l3| . In this situation interlayer in- 
teractions alone can induce correlations between the elec- 
trons in different layers. Moreover, using electric gates 
the doping level in each layer can be controlled indepen- 
dently and can be adjusted to have nesting of the p-type 
Fermi surface in one layer with the n-type Fermi surface 
in the other layer, condition in which the heterostructure 
is said to be balanced. For balanced double layers the 
interlayer interactions can drive the system into an inter- 
layer phase coherent broken symmetry ground state [bH - 
17| . This state can be thought of as an exciton conden- 
sate 3 of electrons in one layer and holes in the other 
layer, as a superfluid state, or by treating the layer de- 
gree of freedom as a spin degree of freedom (pseudospin) 
as a ferromagnetic state. Very recent data [13| provide 
hints that the interlayer phase coherent state might have 
already been realized in symmetric double layer graphene 
systems. In general this state can be realized in a variety 
of setups (see for example [2(J) and there is strong evi- 
dence that it has been observed in quantum Hall bilayers 
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The experimental capability to realize high qual- gates (V„ s , 



ity graphene-hBN heterostructures has opened unprece- 
dented possibilities to study the effects of interlayer elec- 
tronic interactions. In particular it has made possible to 
study the effects of interactions between fermionic quasi- 
particles having both qualitatively different dispersion 
and chirality. This can be realized by creating a het- 
erostructure in which one layer is SLG and the other 
BLG. In SLG the low energy fermionic excitations are 
massless and have Berry phase 7r, whereas in BLG they 
are massive and have Berry phase 2tt 0, ■ In this work 
we study the nature of the interlayer broken symmetry 
state for the SLG-hBN-BLG systems. We find that the 
difference in the dispersion and chirality between the two 
layers profoundly modifies the nature of the ground state. 
In particular we find that due to the difference of chiral- 
ity: (i) the interlayer broken symmetry state is N-fold 
(N=2 or 4 depending on the nature, long-range or short- 
range, of the interlayer interaction) degenerate; (ii) one 
of the degenerate states is always chiral, i.e character- 
ized by a complex order parameter whose phase depends 
on the momentum direction. These facts are unique to 
double-layer systems in which the fermionic excitations 
in the two layers have different chiralities. The N-fold 
degeneracy of the ground state raises the possibility that 
in SLG-BLG system a state could be realized analogous 
to states realized in 3 He (28|. Moreover, the chiral nature 
of one of the ground states makes possible the realization 
of protected midgap states in the presence of vortices in 
the exciton condensate [29l - l3ll |. 

The system that we study is summarized in Fig.[TJ The 
heterostructure is formed by one layer of SLG and one 
layer of BLG separated by a thin dielectric, that in cur- 
rent experimental setup can be assumed to be hBN. The 
two layers (SLG and BLG) are connected to separate 



-V^ as shown in Fig.QJa)) so that their re- 
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Figure 1: (Color online) Schematic of the experimental setup 
(top) and noninteracting bands (bottom), (a) SLG and BLG 
are gated individually at voltages Vg and — Vg 8 , so that the 
total voltage difference between the two layers V g = Vg + V g s 
and the most relevant bands at low voltages are the BLG con- 
duction band and SLG valence band for low-energy physics, 
(b) Inverting the voltages (—Vg and V g s ) so that the most rel- 
evant bands become the SLG conduction band and the BLG 
valence band. 



spective doping can be controlled independently. The 
conduction and valence band in SLG and BLG touch at 
the corners of the BZ. At low energies the band structure 
of SLG is well described by two inequivalent valleys (at 
the K and K' point in the BZ, also called Dirac points, 
DPs)) around which the fermionic dispersion is linear. In 
BLG the low energy conduction and valence bands also 
touch at the points K K' but around these points the dis- 
persion is parabolic with an effective mass m « 0.03 m e . 
We restrict ourselves to the voltage regime in which the 
Fermi energy of BLG is well below the higher energy 
bands (separated by 71 sa 400 meV from the low energy 
bands). In this regime the low energy physics can be ef- 
fectively described by a two-band model formed by the 
conduction band of BLG and the valence band of SLG 
(or vice versa as shown in Fig.[TJb)). For a given bias 
voltage V g = V g B + Vg the Fermi surface in SLG can be 
set to be nested with the one in BLG. 

The low energy physics of the SLG-BLG system is 
described by the Hamiltonian: % = T-Lq + Hint, where 



an electron with momentum k in layer a, c k (ck,<r) 
creates (annihilates) an electron with momentum k in 
layer cr, refers to the intralayer interaction in mo- 
mentum space and V^ the interlayer interaction. As- 
suming, for concreteness, that the gate voltages are such 
that the Fermi energy is in the conduction band for BLG 
(cr = l) and in the valence band for SLG (a = 2) we 



have ek.1 = -V g B + \Jh 2 v 2 F k 2 + \ - y/^ + h 2 v 2 F k 2 ~f 2 
and £k,2 = Vg — hvpk to which correspond the eigen- 



states Vk,i = ^75 (l,e"' 20k 
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and ipk,2 

respectively, where ?y = +1 for states around the K 
point and —1 for states around the K' point. Below 
we consider the states around the K point only as the 
K' point follows similar analysis, vf ~ 10 6 m/s is 
the Fermi velocity of SLG close to the Dirac point and 
6>k = arctan(fc y /fca,). The chiral factors are /i(6*k — P ) = 
\ [l + e 21 ^- "')] , / 2 (0 k -0 p ) = I [1 + e 1 ^- 8 ^] for the 
BLG conduction band and the SLG valence band, respec- 
tively. 

To decouple the interactions we use the mean-field ap- 
proximation and obtain the mean-field Hamiltonian 
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where [A , A = (A x , A y , A z )] are the mean- fields and 



(t x ,t v ,t z )] the 2x2 identity and Pauli ma- 



trices acting in the layer pseudospin space. The trans- 
verse components of the pseudospin effective field Ak 
form a complex order parameter = A£ — iA^, due 
to the interlayer exchange interaction, that characterizes 
the interlayer phase coherent state. The magnitude of 
the transverse field |A^| measures the strength of the 
particle- hole condensate. The mean fields Ao, A are 
given by the following self-consistent equations: 
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with (7 = 1,2 representing the layer degree of freedom 
that can be thought of as a pseudospin and A the area of 
the heterostructure. Ek.a denotes the band energy for 



Ak = 2lE V k -P^( 
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where g = 4 is the total spin and valley degeneracy, e 
the dielectric constant of the embedding media, and d 
the distance between the two layers. The ^^gd term 
is specific to the interlayer Coulomb interaction in the 
direct channel, = l/[exp (£p/feT) + 1] are the 

occupation numbers at temperature T of the renormal- 

ized bands e± = A° ± E k with E k = ^ '(A~) 2 + |A^| 2 , 

and Fi(ftc_ p ), F 2 (0 k _ p ), F d (0 k _ p ) (with k _ p = U - 
6? p ) are angle-dependent chiral-factors arising from the 
wavefunction overlap factors between states V'k.o-, ipp,a-'- 
Specifically, the intralayer chiral factors have the ex- 
pressions Fi(# k _ p ) = j (cos 26* k _ p + cos # k _ p -I- 2) and 
-F2(# k _p) = | (cos20 k _ p — cos6* k _ p ), while the inter- 
layer chiral factor is 



F 



!J i26l k _p^ 



(6) 



Equations ([21)-© differ from the ones valid in the 
symmetric case of double-layer graphene [l4, 15 1 in two 
important aspects: (i) Even when the heterostructure 
is electrically balanced, due to the asymmetry of the 
band dispersion between the two layers the field A£ does 
not vanish; (ii) The chiral-factors -Fi(6> k _ p ), ^(^k-p), 
F d (8\s.- P ) differ from the symmetric case due to the 
difference in chirality between SLG and BLG; in par- 
ticular for the SLG-SLG heterostructure Fg. s (6'k-p) = 
i (2 + e -' te k-p + e ' 0k -p). The latter point causes the 
interlayer phase coherent state of the SLG-BLG het- 
erostructure to profoundly differ from that of SLG-SLG. 

To understand the consequence of the difference in the 
chiral-factor F d (6k- P ) on the gap-equation for Aj^r be- 
tween the symmetric SLG-SLG heterostructure and the 
asymmetric SLG-BLG, let us write the general solution 
of the gap equation (J5J as A^ = |A£ \j e uek+i< >' with the 
chirality J = 0, ±1, ±2, . . . and an arbitrary global phase 
<\>. Without loss of generality we assume A£, A£, and the 
magnitude | A^ | to be angle- independent (it is straight- 
forward to verify that this assumption is consistent with 
the self-consistent mean field equations). From ([5]), we 
obtain the gap equation for the order parameter, 
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Let us first consider the case of short range interactions 
so that = const. In this case we see from (0 that in 
the symmetric SLG-SLG system the solution with J = 0, 
because of the form of the chiral factor Fg_ s (6> k _ p ), has an 
effective interaction that is twice compared to the J ^ 
chiral solutions and therefore the J = ground state has 
a critical temperature higher than that of chiral states. 
On the contrary for the asymmetric SLG-BLG system 
the chiral factor ^ ensures that the effective interaction 



is the same for all the 4 states J = —1, 0, 1, 2. As a con- 
sequence, for the SLG-BLG heterostructure, in the pres- 
ence of short-range interactions the J = — 1, 0, 1, 2 states 
satisfy the same gap equation and therefore, at the mean- 
field level, are degenerate, i.e., the ground state is 4- fold 
degenerate in this case. We can generalize the analysis 
above to the case in which in one layer the electrons have 
chirality m and in the other chirality n (m ^ Q, n =/= 0). 
In this case, for states close to the K point |32| , we have: 



F" 



?k- P J 



gi(m— n)0k-p _|_ e im9 k - 



and from Eq. |J7J| we find the general result that the bro- 
ken symmetry ground state is 4-fold degenerate when 
m 7^ n and non-degenerate when m = n. 

In many cases of interest we can expect that the in- 
teractions will not be short-range but still be central, 
i.e., Vt?_ p = V(|k — p|). For the practically relevant case 

of Coulomb interactions we have V a 



2TTC- 1 



yd _ 

q 



e | q | . In this case the parts on the right hand side 
of Eq. ([7|) that are odd in vanish after integrating 

over angle and the gap equation takes the form: 
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js((n+ J)6» k _p) + cos(j0 k _p) + 
+ cos^(m — n — J)# k - P ^ + cos^(to — J)# k - P ^ 
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From Eq. |[5J we see that for symmetric heterostructures, 
m = n, in the case of "central forces" the J = state 
again has the highest effective pairing strength and there- 
fore the highest critical temperature jlij]- On the other 
hand for asymmetric heterostructures in which m = 2n 
the states J = and J = n (J — and J = — n for 
the other Dirac point) have the same and the strongest 
pairing strength and therefore the ground state is 2-fold 
degenerate. 

For the SLG-BLG heterostructures, in the presence of 
Coulomb interactions, we therefore find that the ground 
state is two- fold degenerate: around the K (K') point, 
the non-chiral J = interlayer phase coherent state 
(layer-ferromagnetic state) is degenerate with the chiral 
state J = 1 (J = —1). This situation is illustrated in 
Fig. [5] where we have chosen the arbitrary global phase 
factor <f> = for convenience. By inverting the gate volt- 
age V g (as illustrated in Fig.QJb)) the values of J at the 
K and K' points are interchanged. 

We find that the nature, chiral or not chiral, of the 
ground state strongly affects the dynamical density- 
density response function for frequencies ui m 2A [33fl 
and therefore that optical measurements should be able 
to distinguish between the two degenerate states. 
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Figure 2: (Color online). Pseudospin configuration on the 
Fermi surface in the broken-symmetry state for a hybrid, 
SLG-BLG, graphene heterostructure around the if-point 
(top) and A"-point (bottom). The J = state, (A£, A£) = 
|A^| (1,0), and the chiral J = 1 (J = -1) state, (A k , A k ) = 
|A^|(cos0 k ,-sin6i k ), ((A£,A fc ) = |A£] (cos6» k , + sin<9 k )) 
are degenerate around the isT-point (if'-point). 



The fact that one of the possible interlayer phase coher- 
ent states is chiral opens the possibility to create topo- 
logically protected midgap states (30l . |34| at the center 
of vortices that can be created in the excitonic conden- 
sate via the axial gauge field, i.e. the part of the vec- 
tor potential A associated to the antisymmetric com- 
ponent Bi — B2 of the magnetic field in the two lay- 
ers [2^|. To see this we observe that we can sepa- 
rate the mean- field Hamiltonian into two parts Hmf = 
Hi + n 2 with Hx = Ek,^ ct a (A k r° CT ,) c k , ff ,, H 2 = 

~J2k.,<r,a' c k,a(^' T ^') c K<r'- Since Hi and H 2 com- 
mute, the eigenvalues of Hmf are given by the sum of 
the eigenvalues of Hi and H 2 . H 2 has a symmetric spec- 
trum ±i?k that in the J = 1 interlayer phase coherent 
state, due to the p x — ip y structure of the order parame- 
ter, in the presence of a vortex in the exciton condensate, 
should guarantee the existence of topologically protected 
mid gap states bounded to the vortex with energy A^ 

We now discuss the conditions for the realization of the 
interlayer phase coherent state in SLG-BLG heterostruc- 
ture. Fig. [3] shows the values of the mean fields as a 
function of |k| (in units of ko = ji/(Hvf)) at zero tem- 
perature, with layer separation d = 1 nm, bias voltage 
V g = O.371, and the effective coupling constant a = 1 
where a = e 2 j {ehvp)- We see that for this parameters' 
set the peak value A = I I of the order parameter 



magnitude is w 0.0757i w 342^. 
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Figure 3: (Color online) Mean fields |A k | , A k , A k , in pan- 
els (a), (b), (c) respectively, as a function of k for T = 0, 
d — lnm, a = 1, and V g = O.371. (d) Renormalized band 
structure (solid lines) of the condensed state. For comparison 
the dashed lines show the noninteracting bands. 

Fig. Ufa) shows the dependence of the order parameter 
amplitude A on the strength of the electric field perpen- 
dicular to the layers for both SLG-BLG and SLG-SLG 
at d = 1 nm and a = 1 . We see that at low electric fields 
(V g /d < 60meV/nm) A is larger in the hybrid SLG-BLG 
heterostructures, whereas at larger electric fields A is 
larger in the symmetric system. Compared to the sym- 
metric SLG-SLG structure, in the BLG-SLG structure 
the density of states (DOS) in one of the layers (BLG) 
is higher than in SLG-SLG structure, and the interlayer 
chiral factor F d (6^ p ) oscillates more rapidly. The first 
effect favors the formation of the interlayer phase coher- 
ent state and therefore enhances A whereas the second 
tends to suppress it. We can then understand the scaling 
of the ratio between A for SLG-BLG and for SLG-SLG as 
a result of the competition in SLG-BLG systems of two 
effects: at low electric fields the DOS effect dominates so 
that A in SLG-BLG is higher than in SLG-SLG, at higher 
electric fields the faster oscillation of F d (8k- P ) causes A 
to be smaller in SLG-BLG than in SLG-SLG. The de- 
pendence of A on the strength of a for both SLG-BLG 
and SLG-SLG at V g = O.271 with d = 1 nm is plotted 
in Fig.[3]Jb). It shows that in the weak coupling regime 
the interlayer coherence becomes stronger in the hybrid 
system at this voltage. 

The value of A for typical gate value of V g w O.371 
suggests a mean-field critical temperature T c > 300 K 
for the establishment of interlayer phase coherence. This 
value is likely an overestimate due to the fact that we 
have neglected screening and disorder effects, and the 
fact that, because the system is two-dimensional and the 
broken symmetry U(l), T c is reduced to the Kosterlitz- 
Thouless temperature Tkt above which we have the pro- 
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Figure 4: (Color online). Dependence of the order- parameter 
amplitude A = |A^ | max on V g , (a), and a, (b), in the hybrid 
SLG-BLG structure and the symmetric SLG-SLG structure 
for T — and d — 1 nm. In (a) a = 1, in (b) V g = O.271. The 
insets show the ratio of A in SLG-BLG to that in SLG-SLG. 



liferation of unbound vortices and antivortices of the con- 
densate. Tkt-, although lower than T c , can be expected 
to be not much lower than T c 



14J. On the other hand 



the screening and disorder effects can in principle reduce 
T c by orders of magnitude. However, it is very difficult 
to correctly quantify the effects of screening and disor- 
der. For example, for the case of the symmetric SLG- 
SLG heterostructures different approximations to treat 
the screenin g eff ects have yielded considerably different 
values of T c \vA [36l - H | . Disorder in these systems also 
induces strong carrier density inhomogeneities [B| whose 
effects on the condensate are difficult to quantify [45[. 
Considering that we find that in SLG-BLG the mean-field 
T c for unscreened Coulomb interaction is of the same or- 
der as in SLG-SLG and that screening and disorder are 
expected to affect T c similarly in the two systems, we 
conclude that in realistic setups T c for SLG-BLG should 
be of the same order as for SLG-SLG. Recent results (l3| 
show hints of an exciton condensate for SLG-SLG in cur- 
rent experimental conditions. We can then conclude that 
the combined effects of screening and disorder in SLG- 
SLG and SLG-BLG heterostructures might suppress T c 
but should not prevent the experimental observation of 
the predicted interlayer phase coherent states. 

The possibility to realize high-quality graphene-based 
heterostructures allows the creation of novel systems of 
great fundamental interest. SLG-BLG heterostructures 
allow the study of the interaction effects among fermionic 
states with qualitative different dispersion, massless and 
massive, and different chirality. We have studied for the 
first time the nature of the spontaneously broken symme- 
try ground state, exciton condensate, driven by interlayer 
interactions in SLG-BLG, and the conditions for its real- 
ization. We find that in SLG-BLG due to the difference 
in chirality between the two layers the ground state is 
2-fold degenerate (4-fold if the interactions can be ap- 
proximated as short range) and that one of the two de- 
generate states is chiral. The chiral nature of the ground 
state should cause the presence of topologically protected 
midgap sates bound to vortices of the excitons conden- 
sate. We find that T c for the formation of the exciton con- 



densate in SLG-BLG is of the same order as in SLG-SLG 
and that at low dopings T c should be higher in SLG-BLG 
than in SLG-SLG heterostructures. The qualitative fea- 
tures of the interlayer phase coherent state of graphene 
hybrid heterostructures, in particular its degenerate na- 
ture, are robust and independent of the details of the 
experimental conditions. 
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